Schwinger Pair Production in Pulsed Electric Fields 
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We numerically investigate the temporal behavior and the structure of longitudinal momentum 
spectrum and the field polarity effect on pair production in pulsed electric fields in scalar quantum 
electrodynamics (QED). Using the evolution operator expressed in terms of the particle and antipar- 
ticle operators, we find the exact quantum states under the influence of electric pulses and measure 
the number of pairs of the Minkowski particle and antiparticle. The number of pairs, depending 
on the configuration of electric pulses, exhibits rich structures in the longitudinal momentum spec- 
trum and undergoes diverse dynamical behaviors at the onset of the interaction but always either 
converges to a momentum-dependent constant or oscillates around a momentum-dependent time 
average after the completion of fields. 
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I. INTRODUCTION 

Schwinger pair production by temporally or spatially 
inhomogeneous electric fields has been a theoretically and 
experimentally interesting and challenging issue. The 
Minkowski or Dirac vacuum under the influence of a 
strong electric field becomes unstable and emits electron- 
positron or charged particle-antiparticle pairs. The pos- 
sibility of directly measuring electron-positron pairs by 
strong laser sources such as Extreme Light Infrastruc- 
ture (ELI) has recently boosted intensive researches on 
pair production by pulsed electric fields [l|-|4j (for review 
and references, see also Refs. 0-0]). In contrast to vac- 
uum polarization, pair production has been studied for 
various electromagnetic configurations, for which many 
analytical approximation schemes have been elaborated 

1- 

The Keldysh approach d, (review and references, 
see Ref. [k|), the worldline instanton method [IJ Ell 
and the phase-integral or WKB method [1343: to list 
a few, have been widely used as analytical approxima- 
tion schemes. The two typical methods employed to 



numerically compute the pair-production rate in time- 
dependent electric fields are the kinetic approach [l9l - [25| 
and the Wigner formalism [26 -28]. The quantum Vlasov 
equation has been often used for numerical calculation of 
pair pr oduction with or without back-reaction included 

The main purpose of this paper is to employ the evo- 
lution operator formalism for numerical computation of 
pair production by pulsed electric fields of various config- 
urations in scalar QED. The stratagem is to express the 
Hamiltonian in terms of the creation and annihilation op- 
erators of the Minkowski vacuum and then find the evolu- 
tion operator satisfying the time-dependent Schrddinger 
equation. In fact, the Hamiltonian for a spinless charged 
boson in a time-dependent electric field in the momen- 
tum space consists of infinite number of oscillators with 
time-dependent frequencies and each oscillator can have 
the evolution operator represented in SU (1,1) algebra. 
We then choose the number operator, the one-pair cre- 
ation operator and the one-pair annihilation operator in 
the Minkowski vacuum as the generators for the evolu- 
tion operator 
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In quantum optics the evolution operator for a time- 
dependent oscillator has been studied in a different repre- 
sentation of SU(1, 1) [37} • The oscillator representation 
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or the creation and annihilation operator representation 
has also been used to derive the quantum Vlasov equa- 
tion for pair production (33M361 HI, [3!| . The difference 
of this paper from other earlier works is that we find the 
evolution operator by directly solving the time-dependent 
Schrodinger equation, then find the exact quantum state 
evolved from the Minkowski vacuum, and finally com- 
pute the number of pairs during the whole evolution in 
pulsed electric fields. The evolution operator represented 
by the creation and annihilation operators has an addi- 
tional advantage of expressing the exact quantum state 
as the squeezed vacuum of multi-pairs of particle and an- 
tiparticle under the influence of the pulsed electric fields. 
Further it provides a good measure for counting the num- 
ber of Minkowski particle and antiparticle pairs. 

To study pair production from nontrivial configura- 
tions of electromagnetic field, we shall consider only 
pulsed electric fields, simplified model fields, which act 
for a finite duration but is uniform over the space. Thus 
the Minkowski vacuum excites into multi-particle and an- 
tiparticle state during the interaction of the electric field 
and finally settles down in some nontrivial vacuum state 
and, in an exceptional case, returns to the Minkowski vac- 
uum. The pulsed electric fields to be studied are classified 
into the mono-polarity type and the di-polarity type. In 
the first class we consider (i) E(t) = E / cosh 2 (t/r), (ii) 
E(t) = E e- t2 / T \ (iii) A,|(t) = (E /u)er t2 ' T \os(ujt), 
and in the second class we also consider (iv) E(t) = 

E /{cOsh 2 ((t-ti)/T)- Eg/ COSh 2 ((t-f2)/ T), (v) A\\ (t) = 

E T/(l+t 2 /r 2 ), (vi) = y/n(n + l)T/ cosh(t/-r). 

The electric fields (i)-(vi) provides a good arena to dis- 
cuss some fundamental questions, not to mention possible 
applications to ultra-strong lasers. The electric fields (i)- 
(iii) of mono-polarity type necessarily lead to non-zero 
gauge potential values, in which the adiabatic basis dif- 
fers from the Minkowski one. Then a question may be 
raised whether pairs produced by electric pulses are the 
Minkowski ones or the adiabatic ones [32l |36|. And the 
oscillation of the number of pairs after the completion of 
electric pulses may be another interesting question since 
the kinetic approach using the adiabatic basis predicts 
non-oscillating pair production while the nonadiabatic 
method predicts oscillating pair production for general 
electric fields except for the solitonic gauge field [36| . The 
evolution operator provides asymptotic solutions in the 
remote future for all pulsed electric fields including (i)- 
(vi) , according to which the number of pairs converges to 
a momentum-dependent constant when the gauge poten- 
tials vanish or oscillates around a momentum-dependent 
time-averaged constant when the gauge potentials have 
nonzero constant. 

Another interesting issue is the structure of the lon- 
gitudinal momentum spectrum of pairs discovered by 
Hebenstreit et al in a sinusoidal field with Gaussian enve- 
lope [2l|. Dumlu and Dunne explained the substructure 
of the spectrum in the inverse square gauge potential by 
the Stokes phenomenon, in which more than two pairs 
of complex turning points for the Fourier mode equation 



contribute either in phase or out of phase to the WKB 
instanton action [l6l . [l7j . The analytical approximation 
schemes such as the Keldysh approach or the WKB or 
phase-integral method without Stokes phenomenon can- 
not explain the substructure of the longitudinal momen- 
tum distribution. The evolution operator formalism is 
efficient enough to compute the longitudinal momentum 
distribution of produced pairs in the pulsed electric fields 
(i)-(vi) and confirms the substructure of the spectrum for 
the critical field strength and Compton time scale. The 
spectra for (i), (ii), (iv) and (v) exhibit rich structures 
but those for (iii) and (vi) have relatively simple struc- 
ture. Bunching of spectrum in the field (iii) is observed 
around one positive and one negative momentum with 
the same magnitude. Further the polarity of the elec- 
tric field (iv) leads to negative momentum dominance or 
positive momentum dominance in the spectrum. 

The organization of this paper is as follows. In Sec. 
II, we introduce the evolution operator in algebraic form 
for scalar QED in a pulsed electric field and express the 
vacuum polarization in terms of the complex parameters 
for the evolution operator. And we express the num- 
ber of pairs by the imaginary part of the parameter for 
the number operator. In Sec. Ill, by solving the time- 
dependent Schrodinger equation, we derive a set of first 
order differential equations for three complex parameters 
for the evolution operator. And we find the asymptotic 
solutions for the pulsed electric fields whose gauge poten- 
tials approach either zero or nonzero constant after the 
completion of the interaction. In Sec. IV we give an intu- 
itive interpretation of pair production by pulsed electric 
fields from the analogy with quantum mechanics. In Sec. 
V, we numerically investigate pair production by several 
configurations of electric fields, some of which have re- 
cently been studied in literature. In Sec. VI, we discuss 
the physical implications of the results and conclude the 
paper. 



II. EVOLUTION OPERATOR APPROACH 

In scalar QED the Fourier-decomposed Hamiltonian 
for a spinless charged boson with mass m in an electric 
field along a fixed direction [36| 



H(t) = 



d 3 k 
(2tt) 3 



4^ +u>fe(t)<^ fc 



Wfc(t) = (k\\ -qA {l (t)) 2 +kl+m 
may have an oscillator representation 
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H k {t) 
where 



2 fi< +) p, + ^(i)(j( + »+j(-)), 



Wfe(i ) 



(1) 



(2) 



(3) 
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and 

N k = a\a k + S_fcSL fc , J { k +) = &\bl k , Jt } = a k b- k - (4) 

Here we have used the oscillator representation in the 
Minkowski vacuum 
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y / 2u> k (t Q ) 



7T k 



— — « - b_fc), 



(5) 



where and b_ k are the particle and antiparticle opera- 
tors. And these operators constitute the SU (1,1) algebra 



[iV fc ,jf>] = ± 2jf\ [jW,J«]=iV fc . 



(0) 



In this paper we shall focus on pulsed electric fields that 
act effectively for a finite duration. Before the onset of a 
pulsed electric field, we may choose a gauge 



A||(t) = 0, E{t)=Q (t<t ), 



(7) 



so that the in-vacuum is the Minkowski vacuum. The 



constant electric field with A 



-Erf will not be con- 



sidered in this paper since no gauge can be chosen to 
make the initial vacuum the Minkowski one, for which 
the in- and out- vacua are defined as asymptotic states. 

The evolution of the Minkowski vacuum follows the 
time-dependent Schrodinger equation 



■ dU k (t) ~ . f M 
i q = H k (t)U k (t). 



(8) 



The evolution operator in time-ordered integral 

U k (t) = T(e~ i/ «o J ** (t ' )dt '), (9) 

does not provide useful information about the quantum 
state unless Q k ~\t) — since the Hamiltonian is then 
entangled such that [H k (t"), H k (t')} ^ for t" ^ t' . The 
particle and antiparticle operators evolve as 



a k (t) = U k (t)a k ul(t) := ^ k {t)a k + v* k tf_ k . 
S-fe(t) - U k (t)b- k Ul(t) := /j fc (t)S_fc + vta\ 

The time-dependent vacuum state is given by 



|0,i)= JJ^(*)|0>in) ) 



(11) 



where |0,in) denotes the Minkowski vacuum. The num- 
ber of the Minkowski particle carried by the time- 
dependent vacuum, which is equal to the number of time- 
dependent particle carried by the Minkowski vacuum, is 

(0,i| a{a k \0,t) = (0,in|o£(t)o fc (t)|0,in) = \v k \ 2 . (12) 



The same is true for the antiparticle production. Further, 
the vacuum persistence is related to the mean number 
401 



|(0,i|0,in)| 2 =exp[-^ln(l + M 2 ) 



(13) 



On the other hand, the SU(1, 1) algebra may lead to 
the evolution operator of the form 



U k {t) 



' e ilk(t)N ke q k (t)J t 



(-) 



(14) 



Here the complex parameters "f k and rj k will be de- 
termined from the evolution equation ©. It should be 
mentioned that Eq. (Q3} differs from Eq. (31) of Ref. [37| 
for a time-dependent oscillator, which would read 



J. 



(+) _ 



Ji 



(0) 



i$l<t>k 
i 



J, 



(-) - .-*t* 



{K k 4> k + 4> k K k + Kl4>\ + <fi\%\). (15) 



Another form of the evolution operator is realized by in- 
verting the Bogoliubov transformation (fTO]) in Ref. [40| . 
The representation ((5J) is particularly useful for pair pro- 
duction since it expresses the out-vacuum as a squeezed 
vacuum of the Minkowski vacuum 



|0,i) = JJe^W^W^'lCin) 

k 

oo 

= n efy * (t) (E^(*)i n *> 



n_fc,in. 



(16) 



where n k and fi- k denote the n particles with momentum 
k and n antiparticles with momentum —k, respectively. 
Thus pair production conserves charge and momentum. 
It further leads to the vacuum polarization 



(0 J t|0 J in)=ex P [-*5^ 7 i(t) 

k 

and the vacuum persistence 

|(0,i|0,in)| 2 =exp[-2^ 7 i fe (0 

k 



(17) 



(18) 



where ji k denotes the imaginary part of 7&. Therefore 
we can find the number of produced pairs from the imag- 
inary part of the parameter j k (t) during the interaction 
of the electric field as well as after the completion of the 
interaction 



\Mt)f 



„27ifc(t) 



1. 



(19) 



III. MASTER EQUATIONS FOR EVOLUTION 
OPERATOR 



The time-dependent Schrodinger equation © together 
with SU(1, 1) algebra leads to the differential equations 
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for the complex parameters 1 



jk 4 


-if] k e- 2 ^ k S,k = 


2 k ' 


(20) 


£,k - 2ijk£k 




2 k ' 


(21) 




Vke- 2 ^ = 


2 fe 


(22) 



The differential equations for the real and imaginary 
parts of complex parameters can be grouped into the set 
that determines pair production 



1 



lik 



-flk 6fc, 



(23) 
(24) 



i 



U = -n{ +) ZRk-\rt-\l + e Rk -tfk); (25) 



and into another set that are relevant for the vacuum 
polarization 



7Rfc 

VRk 



~2 n k +) ~ 2^ ^ Rfc ' 



(-)„-27: 



= -in 



sin(7 Rfe ), 
Ifc cos(7 Rfc ). 



(26) 
(27) 
(28) 



To get the initial data for the parameters, we employ 
the Baker-Campbell-Hausdorff formula to write the evo- 
lution in a single exponential operator 



U k {t) = exp[(l - ilk{t)m(t)J ( k +) + Vk(t)J k ') 
+ (*7fc(*) + \^k(t)vk{t){l - m{t)))N k + ••■], (29) 

where the dots denote polynomials of jk and rjk higher 
than third order. Before the onset of the electric field 
f^l ' = 0, the evolution operator takes the form 



and Q^. (oo) = 2u>k{to) = 2ajfc(oo), the asymptotic solu- 



tions are 



Jik{t) = 7ifc(°°), 

£rAs(*) = |ffc(°o)l cos(2cj fe (oo)i + ^ fc ), 
£ifc(<) = -|Cfc(oo)|sin(2a;A(oo)i + Vfc) 



(32) 



Here three integration constants are tfk and the remain- 
ing two are identified with those at t — oo. This implies 
that the number of pairs per unit volume and per unit 
time approaches constant and the squeezing parameter 
rotates as £ fc = \^ k (oo)\e- l{2u]k{oo)t+ipk \ In the second 
case of more general A||(oo) ^ and f2^~^(oo) 7^ 0, the 
asymptotic solutions are given by 



7ife(*) 

e»(*) 



1 



Ci-^ >(oo) / 



2w fc (oo) yT+cf cos(2w fe (oo)i + z? fe ) 
n£ _) (oo) c 2 + V 1 + c 3 sin(2w fc (oo)t + tf fc ) ' 



^i +) (oo) 



(^M) 2 
(^ +) (oo)) 2 



1-1-2 



(i-fer) 



(4 + V)) S 

1/ 21 



(33) 



where ci, C2 and C3 are integration constants and 
tan#fc = l/c 3 . Further integrating Eq. (|33l) 



1 



7lfc(*) = Cl + 7T ln 



C2 



v/l + c§ sin(2u; fe (oo)f + i? fe ) ,(34) 



we find the number of produced pairs 
N k (t) 



,2ci 



C2 



\Jl + 4 sin(2w fe (oo)f + i? fe ) -1.(35) 



The integration constants of the asymptotic solutions 
(1321) or (J33J) should be determined by solving the evo- 
lution equations (|23|) -()25 |) or may be found from other 
analytical methods. 



Therefore the initial data are given by 



(30) 



£(to) - 7fc(*o) = *?(*o) = 0, 7*(*o) = - 2 n *(*o)- (31) 

For pulsed electric fields such that _E(i) = and 
A11 (t) = constant for large i, we find the asymptotic solu- 
tions. In the first case of An (00) = so that o£ — ^ (00) = 



1 Here we have used the algebraic relation, e XN Me AJV 
E^° =0 M ^ n) i where [N,M] < -°'> = M and [N,M] ( - n '> 

[N, [N, Mjt™" 1 )]. 



IV. SCATTERING PICTURE 

The general feature of pair production can be under- 
stood from the analogy with quantum mechanical scat- 
tering problem. We write the mode equation for the 
scalar field as 

-^-P|fW fc = (m 2 +ki)^, (36) 

where the kinematic momentum is 

P ll (t)=k ll -qA {] (t). (37) 

The mode equation Q36p in the domain of time can be 
interpreted as a nonrelativistic particle with mass 1/2 



and energy e 



kj^ under the potential V(t) 
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— P|(t). We assume that the pulsed electric field has the 
gauge potential Aii(— oo) = but Aii(oo) = constant. 
Depending on the profile of Aii(t) during the interaction 
of the electric field pulse, the potential becomes a barrier 
or well and the asymptotic value V(oo) may be the same 
as V(— oo) or higher (the green line) or lower (the sky 
blue line) (see Fig. 1). The incident wave from the future 
is partially reflected by the well or barrier back into the 
future and partially transmitted into the past 



CtkVk{t) + PkV* k {t) i — > Uk(t), 



where 



Uk(t) 



e —iu k (—oo)t 

^2w/c(-oo) ' 



v k (t) 



e -iui k (oo)t 

^2w fe (oo) 



(38) 



(39) 



First, in the limit of weak field (|«jA.|i| <C to) or large 
momentum such that 



uj k (t) - Wfc(-oo) 



Wfc(-oo) 



«1, 



(40) 



Eq. (|36p is the scattering by a high energy particle over 
a shallow well or barrier. Thus the reflection coefficient 
/3k is suppressed for all ku . Second, in the limit of strong 
field (|gA||| 3> to) or small momentum such that 



u k (t) - Wfc(-oo) 



Wfe(-Oo) 



> 1, 



(41) 



Eq. (|36|) is the scattering over a relatively deep well or 
high barrier. The kinetic energy for small k\\ is smaller 
than for large k» , which means that the reflection coeffi- 
cient for small k» is larger than that for large k\\ . 




FIG. 1: (color online). The schematic plot for the scatter- 
ing problem: the red line denotes the energy of the par- 
ticle, fi 2 — m 2 + k^; the blue line denotes a potential, 
V(t) = — Pi?(t)), where ^4||(=Foo) = 0; the green or sky blue 
line denote a potential, V(t) = —P 2 (t), where oo) = 

but A\\ (oo) ^ 0. 

The above quantum mechanical interpretation may 
have a direct interpretation of pair production when 
A|i(oo) — so that u> k (oo) = uj k {— oo). The positive fre- 
quency solution in the right hand side of Eq. (|3"8"|) splits 



into a positive solution and a negative solution with the 
same frequency in the left hand side. Thus \j3 k \ 2 is the 
mean number of produced pairs by the electric field pulse. 
However, when A||(oo) ^ and uj k (oo) ^= Lo k {— oo), we 
may write the solution in the future by another basis with 
the initial frequency 

(t> k (t) = \i k u k (t) + v k u* k (t) . (42) 

Using the Wronskian 

Wv[u k ,u k ]=i, Wi[v k ,vl]=i, (43) 

we find the relation between two set of coefficients 

fi k = -i(a k Wv[vk,ul] + £ fe WrK,?4])> 

v k = i(a k m[v k ,u k }+ f3 k Wr[v* kl u k ]y (44) 

Thus the number of produced pairs measured with re- 
spect to the Minkowski vacuum 



Wl"-|5(^-/ 



Wfc(-Oo) 

Wfc(oo) 



Oik 




^fc(oo) 

Wfc(-00) 



/ ttffe(-Ooh ^ fcC 2fa,»(oo)t 2 
w fc (oo) / 

(45) 



oscillates with the frequency 2wfc(oo) unless oj k {oo) = 
u>k{— oo), which confirms the the result (1351) from the 
asymptotic solution. Note that |^fc(£)| 2 counted with re- 
spect to the Minkowski vacuum is the same as |/3fc(i)| 2 
counted with respect to the adiabatic vacuum only when 
Wfc(oo) = uj k {— oo). (For discussion on pair production in 
the adiabatic basis in QED, see Ref. [32j.) 

A few remarks are in order. First, the kinetic approach 



0k (t) = a k (t)- 



<j>k(t) = -iwk(t) 



(a k (t)- 



y/2w k (t) 



h{ty 



iS*<»k(t)df , 



y/2w k (t) 



(46) 



applied to the pulsed electric field results in the asymp- 
totic solutions (|3"51) with the coefficients 



a fc (oo) = a k {oo) e - 1 ^ UJ ^ t '^ dt ' +ZUk( - oc '> t \ 
/3 k (oo) = ^(oo)e l -/'- 1 ~" ;fc(t ' )dt '-^ (oo) * 1 . (47) 

where t± is an arbitrary time beyond which io k (t) = 
Wfe(oo). Thus the number of pairs from the kinetic ap- 
proach is the same from Eq. (l3"5|) . that is, \(i k \ 2 = \Pk\ 2 ■ 
Second, some gauge potential may involve an odd func- 
tion An (t) such that for a fixed constant c 



P|l(t) = fc|l +C -A ll (t), 



(48) 
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and Pii (— t) for fey +c is the negative of Py (t) for — (fey + c). 
This implies that the scattering equation for — (fey + c) 
is the backward scattering for (fey + c) in time and has 
the same reflection and transmission coefficients. Thus 
the number of pairs is symmetric in the longitudinal mo- 
mentum around fey = — c in the remote future. 



V. PAIR PRODUCTION IN MONO-POLARITY 
ELECTRIC FIELDS 

For the purpose of numerical calculations, we use the 
Compton unit fi = c = e = m = lso that the time 
is measured in Compton time, the field strength in the 
critical strength and the energy in the rest mass energy 
of the particle: 



tc 



tie 



(49) 



We restrict to the zero-transverse momentum, which is 
equivalent to replacing m 2 = 1 by e\ = m? + = 1 

in all the calculations below. In this unit system n£ 
Eqs. (l23 ]) -([25 ]l read 



2 + PfAt) + k\ 



l + kl 



1 + fej 



(50) 



Further we shall consider only the strong field regime 
since pair production is prominent for an electric field 
near or above the critical strength and dynamical char- 
acteristics is manifest at the Compton time scale. 



E(t) = 



e 

cosh 2 (i) 



The Sauter electric field is a frequently used model 
in QED which is uniform in space but nontrivial in 
time. The Green function and the asymptotic number 
of pairs of the adiabatic particle and antiparticle have 
been known 4l|. The gauge potential is chosen 



A\\(t) = -£ T(l + tanh(- — — )) 



(51) 



such that the initial state is the Minkowski vacuum. The 
duration of the field is effectively characterized by 2t. 
We shift the peak of electric field by 10r only for the 
numerical purpose. 

At the critical strength Eq = 1 for a Compton time 
scale pulse r = 1, the number of pairs is order of one and 
exhibits both the temporal behavior and the substruc- 
ture of the longitudinal momentum spectrum as shown 
in Fig.[5J The pair production for small momentum in the 




FIG. 2: (color online). The number of pairs for the Sauter 
electric field with Eo — 1 and r = 1 is plotted as a function 
of time and longitudinal momentum. 





FIG. 3: (color online). For the Sauter electric field, the num- 
ber of pairs Nk (t) [left panel] and the real and imaginary parts 
of the function £,k(t) [right panel] are plotted as a function of 
time for ku = (red), fen = 1 (blue), fcii = 5 (green), and 
ku = 10 (magenta), where the real part is in the light color 
and the imaginary part is in the dark color. 



Compton unit begins to oscillate around the time of in- 
teraction as shown in detail in Fig.[3j in which the gauge 
potential at t — 10 is the half of the asymptotic value in 
the remote future. The pair production increases soon 
after the electric field is turned on, reaches the maxi- 
mum within a few Compton times and then oscillates 
with large amplitude for small fey and with small am- 
plitude for large fey . The small amplitude oscillation for 
large fey is not shown in Fig. [5] and in the left panel of 
Fig. [3] due to a different order of magnitude. However, 
the asymptotic solutions predict oscillations for all 
fcii since Au (oo) ^ 0, which is numerically confirmed. The 
time-averaged number of pairs monotonically decreases 
as the momentum increases. 

The suppression of pair production for large momen- 
tum is expected from the scattering picture in Sec. IV, 
according to which the particle has a large kinetic en- 
ergy compared to the potential well or barrier and thus 
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FIG. 4: (color online). For the Sauter electric field, the lon- 
gitudinal momentum spectrum of pairs at the time t = 100 
is plotted as a function of momentum in the range of —10 < 
k\\ < 10 [left panel] and magnified in the range of —2 < k\\ <2 
[right panel]. 

has a small reflection coefficient, implying small pair pro- 
duction. The oscillatory and temporal behavior of the 
squeezed vacuum (I16|) can understood from in the 
right panel of Fig-El in which it oscillates with large am- 
plitude for small momentum while it oscillates with small 
amplitude for large momentum, according to the asymp- 
totic solutions (j3"3")l . 

The longitudinal momentum spectrum of pairs shows 
a substructure, which could not be seen by the number 
of adiabatic pairs [4lj . Since the kinematic momentum 
P\\(t) for fci| + 1 and — + 1) is antisymmetric in time, 
the longitudinal momentum spectrum of number of pairs 
is symmetric around fcii = — 1 according to Sec. lIVl which 
is numerically confirmed in Fig.|4j Pairs are minimally 
produced for fc|| = —1, for which P\\(— t) = —P\\(t) and 
LOk{— oo) = ojfc(oo). The pair production is suppressed 
for large momentum as expected. 

t 2 

B. E(t) = E e ^ 
The gauge potential given by the error-function 




leads to the Gaussian electric field 

(t-10x) 2 

E{t)=E e (53) 

The center of the Gaussian field is shifted for the numer- 
ical purpose. The Gaussian electric field decays more 
rapidly and thus produces relatively smaller pairs than 
the Sauter electric field in Sec. IV Al For the numerical 
work we set Eq = 1 and r = 1 as for the Sauter electric 
field. 

The number of pairs as function of time and longitu- 
dinal momentum in Fig. [S] exhibits a similar structure as 




FIG. 5: (color online). The number of pairs in the Gaussian 
electric field with Eo — 1 and r = 1 is plotted as a function 
of time and longitudinal momentum. 




FIG. 6: (color online). For the Gaussian electric field, the 
number of pairs Nk (t) [left panel] and the real and imaginary 
parts of the function (t) [right panel] are plotted as a func- 
tion of time for fcii = (red), fcii = 1 (blue), fcii = 5 (green), 
and fen = 10 (magenta), where the real part is in the light 
color and the imaginary part is in the dark color. 



that of the Sauter field in Fig.[5J It has a structure for 
small momentum but it is suppressed for large momen- 
tum. The similarity of the temporal behavior can be seen 
by comparing the left panel of Fig.[5]with Fig. [31 The lon- 
gitudinal momentum spectrum of pairs in Fig. [7] also has 
a similar pattern as Fig.UJ As the error-function is an 
odd function erf(— t) = — erf(t), it is symmetric around 
k\\ = — y / 7r/2 according to Sec. lIVl which is numerically 
confirmed in Fig. [7] Note that Atm = — v / 7r/2 is the chan- 
nel for minimal pair production for small momentum. 
The suppression of pair production for large momentum 
is also shown in Fig.[7J 
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FIG. 7: (color online). For the Gaussian electric field, the 
longitudinal momentum spectrum of pairs at t — 100 is plot- 
ted in the range of —10 < ku < 10 [left panel] and magnified 
in the range of — 2 < ku < 2 [right panel]. 





FIG. 8: (color online). The profile of electric field (|55ft is 
plotted in red color for Eq = 1, r = 1 and lo — 10 [left panel] 
and lo — 30 [right panel]. The curve in blue color in each 
panel is e - '' -10 ' sin(ajt) for the same value of lo. 




FIG. 9: (color online). The number of pairs in the oscillating 
Gaussian electric field with Eo = 1 and r = 1 is plotted as a 
function of time and longitudinal momentum. 



En, 



s(uot) 



Finally, we consider a gauge potential 

En (t-ioo 2 

An(t) = — e cos(wi), (54) 

UJ 

which leads to the oscillating Gaussian electric field 



E(t) = 



(t-10 T ) z 



sin(cji) 



2Eo 

LOT 



t-Wr 



cos(wi). (55) 



As shown in Fig. [5] the first term is dominant in the re- 
gion \t— 10t\ < 2/ (lot 2 ) and oscillates with an Gaussian 

envelope e~( t ~ 10r - ) l r and beyond that region the sec- 
ond term has a linearly growing factor but both terms 
are exponentially suppressed. 

For the numerical work we set Eq = 1, r = 1 and 
Lo = 10 or Lo = 30. It is surprising that the number of 
pairs in Fig. [S] is bunched around k\\ — —5 and k\\ = 5. 
The longitudinal momentum spectrum in Fig. [10] shows 
bunching with momentum separation of 10 and 30 for 
lo = 10 and lo = 30, respectively. And the rapidly oscil- 
lating factor suppresses pair production, though the peak 



FIG. 10: (color online). For the oscillating Gaussian electric 
field, the longitudinal momentum spectrum of pairs is plotted 
at t — 100 for uj — 10 [left panel] and uj — 30 [right panel]. 



intensity in Fig. [8] is almost comparable to the Sauter 
electric field in Sec. IV Al and the Gaussian electric field 
in Sec. IV Bl The numerical study also shows a similar 
pattern in pair production and bunching of pairs in the 
gauge potential 



En (t-iQ T ) 2 
A\\(t) = — e sin(cjt). 

LO 



(56) 



VI. 



PAIR PRODUCTION IN DI-POLARITY 
ELECTRIC FIELDS 



In this section we study how the polarity of the elec- 
tric field influences pair production, in particular, when 
the electric field changes the polarity during the interac- 
tion. It is interesting to understand how the produced 
pairs behave when another electric pulse of opposite po- 
larity acts. Do they annihilate each other partially or 
completely? Or are there more pairs produced by the 
second pulse? In order to answer some of these questions 
we consider three model fields. 
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E(t) = 



Eg 



Eg 



cos h2(£^i) cosh 2(£^£2) 



As the first model we consider the di-polarity Sauter 
electric field, two Sauter electric fields acting with a time 
lag and in opposite directions. The model gauge field is 
given by 



.4 



-E r (tanh( 



i-lOr, 



tanh( 



t-20r. 



(57) 




Though we cannot exactly solve QED problem in this 
gauge potential, we may understand the characteristic 
feature of pair production by each of Sauter electric fields 
when they are separated by a sufficient time gap as in 
quantum mechanics. 




FIG. 11: (color online). The number of pairs in the di-polarity 
Sauter electric field with Eq = 1 and r = 1 is plotted as a 
function of time and longitudinal momentum. 



FIG. 12: (color online). The number of pairs Nk(t) in the di- 
polarity Sauter electric field is plotted as a function of time 
for feii = 1 (red), fcii = 1.5 (blue), fcii = 5 (green), ku = — 1 
(magenta), ku = —1.5 (gold), and fen = —5 (skyblue) [left 
panel] and the longitudinal momentum spectrum is plotted 
at t — 100 [right panel]. 
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FIG. 13: (color online) 
1 I i 

cosh 2 (t- 10) 1 

fc|| = 1 (red), k 
(magenta), fen = —1.5 



sh a (t 



"II 



The number of pairs Nk(t) in E(t) = 
■jqJ is plotted as a function of time for 
1.5 (blue), fen = 5 (green), fey = — 1 
(gold), and ku = —5 (skyblue) [left 



panel] and the longitudinal momentum spectrum is plotted 
at t = 100 [right panel]. 



For the numerical work we set Eq = 1 and r = 1 
so that two Sauter electric pulses are effectively sepa- 
rated. At the onset of the interaction the peak of num- 
ber of pairs in Fig.[TT] is almost the same as the single 
Sauter electric field in Fig. [2 though it exhibits more 
structure. However, the main difference appears after 
the completion of the interaction. The di-polarity Sauter 
electric field returns to the Minkowski vacuum while the 
Sauter electric field has a constant residual gauge and 
L0k(— oo ) 7^ Wfe(oo). The asymptotic solutions (|52"]l pre- 
dict that the number of pairs for the di-polarity Sauter 
field approaches to a constant, as shown in Figs. [11] and 
[T2l while Eq. (|35|) predicts that the number of pairs 
for the Sauter field oscillates with a constant amplitude 
around the time-averaged value, as shown in Figs. [5] and 

El 

The gauge potential in between two peaks of the di- 
polarity Sauter field has approximately the shape of 
square potential barrier, in which pair production oscil- 
lates according to the asymptotic solution (pi5j) . No pair 
production for k\\ = is the numerical coincidence that 



m = 1 in the Compton unit is equivalent to five wave- 
lengths of a particle in ten Compton length width of the 
square barrier in the scattering picture IIVI and has the 
zero reflection coefficient due to resonance. 

The effect of the polarity of the electric field can be 
seen in the right panel of Fig. [12] and Fig.[T3l in which 
the longitudinal momentum spectrum shows the mir- 
ror symmetry when the polarity of the electric field 
changes. There are residual pairs with negative momen- 
tum (Fig. [12]) when the first Sauter field acts in the pos- 
itive direction and then the second one in the negative 
direction, while more pairs with positive momentum sur- 
vive (Fig. [T3]) when the first Sauter field acts in the neg- 
ative direction and then the second one in the positive 
direction. The number of pairs in time also shows the 
polarity effect as shown in the left panel of Fig.[T2] and 

Fig.E!J 
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FIG. 14: (color online). The number of pairs in the inverse 
square potential is plotted as a function of time and longitu- 
dinal momentum. 





FIG. 15: (color online). For the inverse square potential, the 
number of pairs Nt (t) [left panel] and the real and imaginary 
part of the function £k (t) [right panel] are plotted as a function 
of time for ku — (red), ku = 1 (blue), k» = 2 (green), and 
fcii = 5 (magenta). 



The inverse square potential shares the same property 
of vanishing in the remote past and future and chang- 
ing polarity as the di-polarity Sauter electric field in 
Sec. lVI Al The structure of the longitudinal momentum 
spectrum of pairs in the inverse square potential was dis- 
covered and explained as a consequence of Stokes phe- 
nomenon by Dumlu and Dunne [lg, Qj} ■ 

For the numerical work we set Eq = 1 and r = 1. 
The peak of the number of pairs in Fig.Q3] is relatively 
small compared with the Sauter or Gaussian electric 
field almost with the same peak intensity. Also it is 
small even compared with the di-polarity Sauter electric 
field. After the completion of interaction the number of 
pairs oscillates with small amplitude as shown in Fig. [15] 
The structure of the longitudinal momentum spectrum is 
shown in Fig. 1161 The relatively simple looking spectrum 
has also a fine substructure as shown in the right panel 



of Fig.|T5J The result of this section cannot be directly 
compared with that in Ref . [l6l Il7j , in which the struc- 
ture was found for a subcritical strength and longer time 
scale than ours in the Compton scale. 



C. Adt)= V^^ L 

II v ' cosh( — J 



The solitonic gauge potential has a special energy con- 
dition 

, . . n, n n(n + 1)cj 2 (— oo) 
u J 2 (t)=ui 2 (-oo) + K > K ' (58) 
cosh oojij 

where n is a natural number and oj(— oo) is the Minkowski 
energy of particle in the remote past. The one-soliton 
(n = 1) and multi-soliton (n = 2, • • • ) gauge fields corre- 
spond to soliton solutions for the inverse scattering of the 
Korteweg-de Vries (KdV) equation [3(|[42j]. The solitonic 
gauge potential is interesting in understanding the pecu- 
liarity of Schwinger mechanism by pulsed electric fields. 
Among the gauge potentials that vanish in the remote 
past and future, the solitonic gauge potential is unique 
in that the number of pairs produced by the correspond- 
ing electric field exponentially increases from and then 
decreases to zero [H, H3|. In fact, any gauge potential 
which has the reflectionless scattering for the field equa- 
tion (|36p has zero number of pairs in the remote future 
and may belong to the solitonic gauge potential. 

For the numerical purpose, we consider the electric 
field with the center shifted 



E(t) 



cosh 2 



(59) 



and set r = 1 but select a natural number and non- 
natural number for the strength Eq. The left panel of 
Fig. [IT] is the number of pairs as a function of time for 
the solitonic gauge potentials Eq = 1,2,5 and 10. The 
number of pairs increases and then decreases in symmet- 
ric way around the center t = 10 of the gauge potential 
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FIG. 17: (color online). The number of pairs N(i) is plot- 
ted for for the solitonic gauge potential with Eq = 1,2,5, 10 
[left panel] and for non-solitonic gauge potential with Ea = 
0.2,0.5, 1.5 [right panel]. 



and the number of local maxima is the same as the soli- 
ton number [42j. On the other hand, for a non- natural 
number, for instance, Eq — 0.2,0.5 and 1.5, the number 
of produced pairs increases from zero but decreases to a 
constant value. This temporal behavior of the number of 
pairs is expected from the asymptotic solutions (|32|) for 
u){— oo) = w(oo) regardless of Eq. It should be recollected 
that the inverse square potential has also finite number 
of pairs in the remote future in Sec. lVIBl in contrast to 
the solitonic gauge field. 



VII. DISCUSSION AND CONCLUSION 

We have employed the evolution operator formalism 
to numerically study pair production in scalar QED by 
pulsed electric fields. The pulsed electric field is nontriv- 
ial in that it acts for a finite period of time and has the 
inhomogeneity of time. The Hamiltonian for a spinless 
charged boson in time-dependent electric field is an infi- 
nite sum of oscillators with time-dependent frequencies. 
After expressing the oscillator Hamiltonian in terms of 
the creation and annihilation operators of particle and 
antiparticle in the Minkowski vacuum, the evolution op- 
erator (fT4|) in SU(1, 1) algebra is factorized into the pair 
annihilation part, the number part and the pair creation 
part. The advantage of this factorization is that the exact 
quantum state (|16p from the Minkowski vacuum under 
the influence of a pulsed electric field is the squeezed vac- 
uum of particle and antiparticle with a complex phase 
factor. The sum of all complex time-dependent phase 
factors for the number operator determines the vacuum 
polarization and the vacuum persistence, which in turn 
is related to the number of pairs produced by the pulsed 
electric field. 

Now the time-dependent Schrodinger equation is 
equivalent to a set of first order differential equations for 
the complex parameters for the evolution operator. The 
evolution of the Minkowski vacuum is governed by three 
parameters in Eqs. (j23]) -(|25 p . and thus the set of first or- 
der differential equations can implement numerical works 



for pair production in various configurations of the elec- 
tric fields. Remarkably the differential equations have the 
asymptotic solutions ([32]) when the gauge potential van- 
ishes and another solutions (pJ5|) when the gauge potential 
approaches a constant value. These asymptotic solutions 
put a strong constraint on the behavior of number of pairs 
such that the number of pairs per unit volume and per 
unit time is a constant for the zero-gauge potential in the 
remote future but it oscillates around the time-averaged 
value for a non-zero gauge potential. 

For numerical works we have selected two classes of 
electric fields or gauge potentials. In the first class the 
electric field does not change the polarity of the field. 
Therefore, the gauge potential, as the negative of the 
time integral of electric field, should have a non-zero 
value in the remote future. We have considered (i) Sauter 
electric field, (ii) Gaussian electric field and (iii) oscillat- 
ing Gaussian gauge potential as the first class of mono- 
polarity electric field. In the second class the electric field 
changes the polarity and has the gauge potential which 
vanishes in the remote future. We have considered in 
the second class (iv) di-polarity Sauter electric field, two 
Sauter fields separated by a time gap, (v) inverse square 
gauge potential and (vi) solitonic gauge potential, all of 
which have the di-polarity of field. 

We have computed the number of pairs as a function of 
time and longitudinal momentum for the pulsed electric 
fields of (i)-(vi) with the Compton scale. In the first 
class (i)-(iii), the number of pairs increases when the 
field acts and then oscillates around the time-averaged 
value as predicted by the asymptotic solutions. In gen- 
eral the number of pairs is suppressed for large longitu- 
dinal momentum as expected from the field equation in 
Sec. lIVI The structure of longitudinal momentum spec- 
trum is found for small momentum for (i)-(iii). Heben- 
streit et al found the substructure of the spectrum in a 
sinusoidal electric field with Gaussian envelope [21] while 
the sinusoidal gauge potential is considered in this paper. 
Also, the structure of the momentum spectrum has been 
found for the class (v) by Dumlu and Dunne [TH, [l7{ . 
The number of pairs as a function time or a function of 
momentum shows similarity between the Sauter electric 
field and the Gaussian electric field. However, the mo- 
mentum spectrum for the oscillating Gaussian electric 
field reorganizes and bunches around one positive and 
one negative momenta with the same magnitude and the 
separation in the momentum space equals to the angular 
frequency of the oscillating electric field. We do not have 
a simple physical explanation for this bunching effect due 
to the oscillating field with a Gaussian envelope. 

In the second class of di-polarity field (iv)-(vi) in which 
the gauge potential vanishes in the remote future, the 
number of pairs increases and then decreases to a con- 
stant rate and has relatively simpler structure of the lon- 
gitudinal momentum than mono-polarity electric fields. 
Still the longitudinal momentum spectrum for the di- 
polarity Sauter electric field exhibits a surprising fea- 
ture of negative or positive momentum dominance of 
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produced pairs depending on whether the positively di- 
rected pulse or the negatively directed pulse acts first, 
which is then followed by the oppositely directed sec- 
ond pulse. Further the solitonic gauge potential produces 
pairs in a symmetric way in time and returns back to the 
Minkowski vacuum without any residual pairs. 

In this paper we have confined our study to scalar QED 
in pulsed electric fields. However, there are a few issues 
related to but not treated in this paper. First, the evo- 
lution operator formalism can also apply to spinor QED, 
since the spin-1/2 fermions have SU(2) algebra isomor- 
phic to SU(1, 1) and may have an evolution operator sim- 
ilar to the spinless boson. Another issue is the vacuum 
polarization, the real part of the in- and the out-vacua 
scattering amplitude. To get a finite effective action one 
should properly regularize the sum over the momentum 
of the real part of the scattering amplitude. Still an- 
other issue is the back reaction of produced pairs. The 
strong field investigated in this paper produces roughly 
one pair per unit Compton volume, which could gener- 
ate an induced electric field screening the external field 
and cause another mechanism for the oscillation of pro- 
duced pairs [1^, |3(| • Finally, it would be extremely useful 
and promising for ELI experimentation to have a simi- 



lar evolution operator formalism for both spatially and 
temporally localized electric fields. All these issues go 
beyond the scope of this paper and will be addressed in 
future publication. 
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